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D ■ Abstract 



Mechanism of self-organization in unbounded, double-species, two-dimensional (2D) point vor- 
tex system is discussed. A kinetic equation is obtained using the Klimontovich formalism. No 
axisymmetric flow is assumed. The obtained collision term consists of a diffusion and a drift term 
similar to the Fokker-Planck type collision term. It is revealed that a mechanism for the 2D inverse 



T3 . 

cascade is mainly provided by the drift term, as the sign of the drift term for the negative vortices 

o 

is opposite from the one for the positive vortices. When the system reaches a quasi-stationary state 

near the thermal equilibrium with negative absolute temperature, dw/dip is expected to be posi- 

> ■ 
\r, tive, where u is the vorticity and t/j the stream function. In this case, the diffusion term dissipate 

O 

^T) ' the mean field energy, while the drift term produces the mean field energy. As a whole, the mean 

^; ■ field energy is conserved. Similarly, the diffusion term increases the entropy, while the drift term 

*y-s decreases the entropy. As a whole, the entropy production is positive (H theorem). It ensures that 

the system relaxes to the global thermal equilibrium state. 



PACS numbers: 47.32. C-, 47.27.tb, 05.10.Gg, 05.20.Dd, 52.27.Jt, 47.27.-i 



I. INTRODUCTION 

The two-dimensional (2D) point vortex system has been successfully applied to under- 
stand the various phenomena including 2D turbulence ll 3|, neutral |4| and nonneutral 

nn 

[5|, 16[ plasmas. These phenomena share a common keyword, "inverse cascade" . To tackle 
the problem, Onsager introduced a concept of negative temperature for the 2D point vortex 
system [7|. To understand the equilibrium state of the 2D point vortex system, or 2D flow, 
much research effort has been devoted. Joyce and Montgomery presented the sinh-Poisson 
equation which gives a thermal equilibrium distribution for the inviscid Euler equation [8(. 
However, surprisingly a time-asymptotic distribution traced by the decaying Navier-Stokes 
equation can be described by the sinh-Poisson equation [9( . Thus, another essential approach 
has been made to understand the relaxation process toward the thermal equilibrium state 
through the kinetic theory. 

Collision term for the 2D single- and multi species point vortex system was presented, for 



multi specK 
?il [a, LoUli 



example, by Chavanis, Dubin and O'Neil |5|, llOH14|. However, these results have a problem 
that an axisymmetric, monotonically decreasing vorticity profile and an angular velocity 
profile do not relax towards the Boltzmann distribution. 

In this paper, we discuss a mechanism of the self-organization in unbounded, double- 
species 2D point vortex system. We extend the kinetic theory of Yatsuyanagi [15( based on 
the Klimontovich formalism. A kinetic equation is obtained which has the Fokker-Planck 
type collision term consisting of a diffusion and a drift term. To derive an explicit formula 
of the collision term, the Klimontovich formalism is used. The obtained collision term 
has following good properties that resolve the above issues; (a) it conserves the mean field 
energy, (b) it satisfies the H theorem, (c) it vanishes when the system reaches a global 
thermal equilibrium state. The H theorem ensures that the system relaxes to the global 
thermal equilibrium state. It is also revealed that the drift term plays an important role in 
the 2D inverse cascade process. When the system reaches a quasi-stationary state near the 
thermal equilibrium with negative absolute temperature, du/dip is expected to be positive, 
where u is the vorticity and ip the stream function. In this case, the diffusion term dissipate 
the mean field energy, while the drift term produces the mean field energy. Similarly, the 
diffusion term increases the entropy, while the drift term decreases the entropy. 

The organization of this paper is as follows. In Sec. II, the point vortex system is briefly 



introduced. In Sec. Ill, we demonstrate explicit formulae for the diffusion and the drift 
terms as intermediate results. In Sec. IV, a detailed calculation of the diffusion term is 
shown. As the similar calculation can be applied to the drift term, details for the drift 
term are omitted. In Sec. V, three good properties of the collision term are demonstrated. 
Finally in Sec. VI, we discuss the important role of the drift term in the 2D inverse cascade 
with negative absolute temperature. 

II. POINT VORTEX SYSTEM 

Let us consider a 2D system consisting of N + positive and N_ negative point vortices. 
The circulation of each point vortex is given by either Q or —Q where Q is a positive constant. 

Cj = Cj + + u>_, (1) 

N+ 

u + = "£n5(r-ri), (2) 

N + +N- 

u- = - £ QS(r-n). (3) 

i=N + + l 

The position vector and the strength of the ith point vortex is given by r, = Tj(t) and 
Qj, respectively, Cj = Cj(r,t) is the ^-component of the microscopic vorticity on the x — y 
plane, and S(r) is the Dirac delta function in two dimensions. For brevity, we shall omit 
the dependences on t and r, if there is no ambiguity. The discretized vorticity (T2]) and (J3]) 
are formal solutions of the microscopic Euler equations. The microscopic variables in the 
microscopic equation are identified by : . 

— uj+(r, t) + V • (u(r, t)u+{r, £)) = 0, (4) 

— d;_(r, t) + V • (tt(r, t)u.(r, t)) = 0. (5) 

Note that the motion of the positive vortices is affected by the negative one through the 
velocity ii appearing in the convective term in Eq. (J3J), and vice versa. From now on, we 
shall denote the two equations for the positive and the negative vortices into the single 
formula with double-sign like: 

— u±(r, t) + V • (u(r, t)u±(r, t)) = 0. (6) 



Other microscopic variables are defined by 

u = u(r,t) = — z x V?/>, (7) 

$ = -0(r,t) 
= fdr'G(r-r')w{r',t) 

= EW-n), (8) 

i 

G(r) = ~ln|r|, (9) 

Z7T 

where it and ■?/> are the velocity field and the stream function in the 2D plane, z is the unit 
vector in the ^-direction, and G(r) is the 2D Green function for the Laplacian operator with 
an infinite domain. 

The microscopic Euler equations ((6]) are formally identical with the macroscopic Euler 
equation 

-u{r, t) + V ■ («(r, t)cu(r, £)) = (10) 

where u(r,t) = (u(r,t)) and u(r,t) = (u(r,t)) are the macroscopic vorticity and the 
macroscopic velocity, respectively. Operator (•) is an averaging operator. As a solution of a 
macroscopic fluid equation should be given by a smooth function, the singular solutions ([2]) 
and ([3]) should be regarded as not the solution of the macroscopic equation but the one of 
the microscopic equation. Thus we call the equation that has the microscopic point vortex 
solution, "microscopic" Euler equation. 

Startin g eq uation is the microscopic Euler equations ([6]). Following the Klimontovich 
formalism 16( , the microscopic variables are decomposed into the macroscopic part and the 
fluctuation. 

Cj± = u± + 5u±, (11) 

u = u + 5u, (12) 

Inserting Eqs. ffTTl) and (fT2l into Eqs. (jSJ), and taking the ensemble average, Eqs. (jSJ) are 
rewritten as the following macroscopic equation with the diffusion fluxes T± = T±(r,t) 

^u ± + V • («w±) = -V • r±(r, t), (13) 

T± = (SuSu±) 

= - J ' dr'F{r - r')(5u'5u±), (14) 

F(r) = zx VG(r). (15) 



We note Su' for 8u(r',t). Similarly, we shall note u/ for u(r',t). The following relation is 
utilized. 

5u = - f dr'F(r - r')5u' (16) 

In the next section, we will analytically assess the diffusion fluxes T±. 

III. EVALUATION OF DIFFUSION FLUX 

To evaluate the diffusion fluxes T± explicitly, we introduce a small parameter e. We 
consider a point vortex system with large N± keeping the total circulation N±Q constant. 
The magnitude of e is of the order of either 1/N± or Q. Orders of the other quantities are: 



± 



0(e> 



du du 

_^0 (e) , -«0(e) 

u « V 2 </> « 0(e°), Vu « 0(e), 

w « V^ w 0(e°), Vix w V 2 ^ w 0(e°), 

Wit w 0(e). 



The expansion parameter e is the same one as Chavanis in Refs. 10|, ll3] and the references 
therein. Expressing T± in the form of a perturbation expansion and gathering the second 
order terms, analytical formulae for the diffusion fluxes r± will be obtained. 

To rewrite the diffusion fluxes in Eq. ffl3|) according to the above prospect, we introduce 
linearized equations obtained by inserting Eqs. (TTTT) and ( I12J1 into Eqs. ([6]) and assembling 
the first-order fluctuation terms: 

d _ 

— 6u± + V • (u6u±) = -8u ■ Vw±. (17) 

at 

As the macroscopic quantities u appearing in the second term in the left-hand side and 
Vw in the right-hand side are supposed to be constant in the time scale of the microscopic 
fluctuation, Eqs. (TT7|) can be integrated: 

5u± = — dr5u(r — u(t — t),t) ■ Vtu± 

Jto 

+6u±(r-u(t-t ),t Q ), (18) 

5u±=- [ dr8u(r' - u'(t -t),t)- Vu/± 

Jto 

+8co±{r'-u'{t-t ),t ), (19) 



where V'u;± = V r .'o;±(r / ). We call this approximation "straight-line approximation" where 



the trajectory of the point vortex is straight 12]. The value of to is chosen to satisfy 
t — t 3> t c where t c is a correlation time of the fluctuation. Substituting Eqs. (ITS]) and ( TT§j) 
into the correlation terms in Eqs. (TH|) . we obtain 



/ (- C drbu[r' - u'(t -t),t)- Vu') 

x (— / dr5u(r — u(t — t),t) • Vu;± 

+ ((- [ dr5u(r' - u'{t -t),t) ■ V'w'] 6u±(r - u(t - t ),t ) 

+ ( 5w(r' — u'(t — t ), t ) ( — / dr5u(r — u(t — t),t) ■ Vtu± 

+ (Su(r' - u'(t - to),to)8u±(r - u{t - t ),t )) (20) 

— / dr (Su(r — u(t — r), t)5u>') ■ Vu± 
J t 

- T dr (5u{r' - u\t - r), t)5uj±) ■ VV + 0(e 3 ) (21) 
J to 

I dr I dr"F[r - uit - r) - r") ■ Vw± (5u{r", r)Su') 
+ j dr j dr"F(r' - u'{t - r) - r") ■ W (^(r", r)5w ± ) . (22) 



When obtaining formulae (1211) . we assume that the first terms in formulae (1217]) are negligible 
as they have two nablas. We drop the last terms also as they should have a factor of l/(t— 1 ) 
and we focus on t — to ^> t c case. The time is shifted from to to t using the straight-line 
approximation. When rewriting formulae (I2T]) as (|22j) . Eq. (Tl6l) is used. 

Inserting Eq. (1221 into Eqs. (TH|) . the following intermediate results are obtained: 



r± 



- f dr f dr' f dr"F{r - r')F(r - u(t - r) - r") ■ Vw± 
Jt J J 

x (6u(r", t)Su/) 

- I dr I dr 1 j dr"F{r - r')F(r' - u\t - r) - r") ■ VV 
Jto J J 

x(5u(r",T)5u ± }. (23) 



IV. EVALUATION OF THE FIRST TERM OF THE DIFFUSION FLUX 

There are three correlation terms (Su)(r",r)Su'), (5u(r",T)5u + ) and (5u(r",T)5u_) in 
Eqs. (|2"5|) . As the expression of the correlation terms (Scu(r",r)Scu + ) and (Su)(r",r)Su)J) are 
very similar to (5u(r",r)5u>'), only the detailed derivation for (5cj(r",T)5u)') will be shown: 

(6u(r", t)8lo'} 

= ((u>+(r", r) + w_(r", r) - co + (r", r) - w _(r", r)) 

x [Cj' + + Cj'_ — u' + — u'_\ \ 
= ((Cu + (r",T)+u4r",T))x(uj' + + &_)) 

-(u + (r",T)+co.(r",T))x(co' + + co'_) 

lN + +N^ \ 

= ( E ^(r" - r l (r))5(r' - r.^ 
+ ( E E nin^(r"-r i (r))(J(r , -r i )) 

\ i=l j^i I 

-u(r",r)u'. (24) 



The first term in the last result in Eq. (1241) corresponds to the case of i = j, and the second 
term corresponds to the case of i ^ j. 

For the i = j case, the formula is rewritten as 

lN++N^ \ 

( E tf5(r"-r t (T))5(r'-r t )j 

/N++N- v 

= / E ti$8{r" - n(r) - r' + ri )5{r' - r,)\ 

7V++7V_ 

= E ^(^"-n^-r' + r^r'-n)). (25) 

8=1 

Here we introduce a stochastic process to evaluate r, — t*j(t): 

Vi - n(r) = J u(r(r'), r')dr' + £ 

« u'(t - r) + £. (26) 

The first term in Eq. (1261) represents the straight-line approximation and the second term 
represents a Brownian motion. The stochastic process represented by (-)^ includes all the 



possible motion to reach position r, at time t. Then, Eq. (]25l) can be rewritten as 

N++N- 

J2 S% (6(r" - n(r) - r' + n)6{r> - Pi )> 

1=1 

Af + +Af_ 

= £ ^<5(r''-r' + W '(t-<r)H-£)}<5(r'-n)> 

8=1 

= (<J(r" - r' + u'(t - r) + £)> € ^K " <"-) (27) 

For the i ^ j case, we introduce an approximation that the correlation between the 
particles can be neglected 

N++N- N++N- 

E E ^({(r'-r^Mr'-rj)) 

AT++JV- JV++JV- 

« E E «i«i <<V - n(r))} <5(r' - r,-)> . (28) 

i j^i 

Also we assume the followings: 

N++N- 

E n t (5(r"-r z (T)))=tu + (r",T)+tU-(r",T) (29) 

i 

u+{r' , ,T) = N + n(8{r"-r i {T))), (30) 

w _(r",r) = -iV_^(5(r /, -r i (r))). (31) 

Inserting Eqs. ([30]) and (JSH) into Eq. ([28]), we obtain 

AT + +7V_ N++N- 

E E 0,0,- <5(r" - n(r))> <5(r' - r,)) 

/ iV++A r _ \ /iV++Af_ 

£ a t <5(r" - r< (r))) x £ fy (5(r' - r,)> 

'N+ \ / 7V+ 

Ea(5(r"-r,(r))) x £ ^(^-r^ 

(N+ \ /N++N- 

+ E «i (^" - r«(T))> x £ fy (5(r' - r,)> 

/N++N- \ /JV+ 

+ E «* W»" ff " n(r))> x £ O, (5(r' - r,)) 

\i=JV + +l / \j=l 

/N++N- \ / N + +N- 

+ E «* (<*(**" - ^(r))> x E «i <*(*•' - *v)> 

\i=7V + + l / \j=N++l& 



iv + o -" ' 7 (iv + - i)n 



N.tt y T ' JV+n 



+N_n u - (r "> T \N- - m^- 

N_tt v ; N.tt 

1 u+(r", t)u>' + - -U;_(r", r)u'_. (32) 



N + TV ' ' + AT. 
Combining the results of the £ = j and the z 7^ j cases, we rewrite Eq. (|24p as 

(<Jw(r", r)<Jw') 



fi <<J(r" - J" 7 + u'(* - r) + £)) f (o^ 



u/ 



^ + K,^;-^>",r)^). (33) 



Similarly, we obtain 



(8u(r",T)Su ± ) 

= ±n (S(r" -r + u(t-r) + £))$ w ± 

-^±(rV)«±. (34) 

At first, we focus on the evaluation of Eq. ( I33p . After that, we will return to the 
evaluation of Eqs. ([3"4"]) . As we request the total circulation N±Q constant, the three terms 
in the right hand side of Eq. ( 133]) are of the same order. To proceed with the evaluation of 
Eq. ( 133]) . a conservation law is introduced 

fdr'(5uj{r",T)5u') = 0. (35) 

Inserting Eq. ( 13"3"]) into Eq. (1351) . we obtain 

j dr'Vt {5(r" -r' + u'(t - r) + £)> s « " w -) 

1 -w+(r", r)wl - ^w_(r", rV 



AT, TV y + M 



w' 



n I dr' (6(r" -r' + u'(t-r) + £)) s (u\ 
-^-c + (r", T )fu,' + dr> - ^_(r", r)Jj_dr' 
Q j dr' (5{r" - r' + u\t - r) + g)) ( (u' + - u/ 
1 -u + (r", t)N + Q - — w_(r", t)N^Q 



0. (36) 



This equation yields 

w±(r",T) 
= ±Jdq' (6(r" -q' + u(q')(t - r) + £)) { co ± (q', t) (37) 

where dr' is replaced by dq' to avoid ambiguity. This equation enables that all the quantities 
at r is converted by the ones at t. Inserting Eqs. (133]) and (j3"7|) into the first term in Eq. 
( 123|) . we obtain 

- J dr j dr 1 f dr"F(r - r')F(r - u(t - r) - r") • Vw± 
x (<Jw(r", r)5a/> 
= -Q [ dr [ dr' [ dr"F(r - r')F(r - u(t - r) - r") • Vw± 
x (5(r" -r' + u'(t - r) + £)>£ (u' + - w'_ 



+ f dr f dr' f dr"F(r - r')F(r - u(t - r) - r") • Vw ± 

x ^«; | <V (5(r" - g' + u(q')(t - r) + £)> s w+C^, i) 

-^'_Jdq'{5{r"-q' + u{q'){t-r) + i))^{q\t) 
We proceed with the evaluation of the second term in Eqs. ( 1381) : 
/" dr / dr' I dr"F(r - r')F{r - u{t - r) - r") • Vw± 

1*/+ I «# (8(r" -q' + u(q')(t - r) + ^ u + (q J , t) 

- ^'_jdq> (S(r" -q' + u(q')(t - r) + Z))^_{q\t) 

j dr'F{r - r') f dr f dq' (F(r - q' - (u - u(q'))(t - r) + £)} ? ■ Vw± 



(38) 



x 



— w + (g / ,t)w+ - —u-(q , ,t)u'_ 



(39) 



o?fc z x ik „ 

— -Vu; ± 



= fdr'F(r-r') [ dr f , , 
J y J Jt J (2tt) 2 |fe| 

x / dq' exp [ik ■ {r — q' — {u — u(q')){t — t))\ (exp(ifc • $ t )) i 

To rewrite formulae (139]) as (HOj) . we use the Fourier transformation: 
F{r-q'-{u-u{q')){t-r)) 



(40) 



10 



12*)* J dk^^eMik-(r-q'-(u-u(q'))(t-T))). (41) 

The term (exp(ik ■ £)), represents a Brownian motion of the point vortices with diffusion 
tensor D and is evaluated by the cumulant expansion: 



(exp(ik ■ £)) € = exp f (i - r) J 

= exp(— v(t — r)) 

where f is a small positive parameter. Inserting the following formula into Eqs. ( l4"Uj) . 

/ <irexp[— ?fc ■ (it — u(q') — v)(t — r)] 

« 7T(5(fc ■ (it — it(q'))) 

we obtain 



(42) 



ik ■ (it — u(q')) 
\k ■ (u - u(q'))\ 2 + v 2 ' 



(43) 



, ., , , i , / Qjfxi Z y\ IKi m , 



x / dq' exp [ifc • (r — q' — (it — u{q')){t — r))] (exp(ifc • £)) £ 



—u+iq^t)^ - —u_tf,t)u'_ 



fdr'F{r-r') f 
x /rfq' 



rffc z x ik 
(2vr) 2 |fc| 2 



<J(fe • (u - u( 9 ')) 



ik ■ (u — u(q')) 

\k ■ (u - u(q'))\ 2 + v 2 



exp(ik ■ (r — q')) 



x 



jj-u+{q', *)w+ - jj-w-{tf,t)u'_ 



(44) 



We substitute r + q" for q' and expand it(q') and u)(q') in the form of the Taylor series 
and retain the zero-th order terms only: 



u(q')=u(r) + q"-Vu(r) + 0(e), 
u(q') = u(r + q") = iu(r) + 0(e). 

Inserting Eqs. (1451) and (1461) into Eqs. (1441) . we finally obtain 
r , . . .. r dk z x ik ^ 



(45) 
(46) 



x /dg' 



(27 

7T(5(fc ■ (it — u(q')) 



ik ■ (u — u(q')) 
\k ■ (u - u{q'))\ 2 + v 2 



exp(ifc ■ (r — q')) 



—u + (q',t)uj' + - —w-{tf,t)w'_ 
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f dr'F{r-r') f 



dk z x ik 

(2vr) 2 |fc| 2 



x / dq" 



7i5(-k ■ (q" ■ V)it) 

1 



|ife- (q" • VH 2 + z/ 2 



exp(— ik ■ q") 



x r^ 



iV_ 



-CU_(X? 



(47) 



It is found that by substituting A; = —A; and q" = —q", the sign of Eqs. ( 147|) changes. 
Thus it is concluded that the integral equals zero, i.e. the second term in Eq. ( 138]) has zero 
contribution and only the first term remains. The obtained formulae for the first terms of 
the diffusion fluxes ( |23|) are as follows: 

- f dr fdr' f dr"F(r - r')F(r - u(t - r) - r") ■ Vw± 

x (5u(r", t)5u'} 



n I dr'F(r - r') f -^ exp(ik ■ (r - r')) 



, z x ik 



\k\ 



w. 



wM Vw± 



7T(5(fc • (it — It')) 



ik ■ (u — u') 



(48) 



|fe • (it — it') | 2 + v 2 

The similar calculations can be adapted for the remaining two terms with (5u(r", t)5uj + ) 
and (Soj(r",r)5oj-). For these cases, the following conservation laws are used: 



dr(Su(r",r)5w±(r,t)) = 0. 



The whole results including both the diffusion and the drift terms are given by 



(49) 



X 



7T(5(fc ■ (it — It')) 



ifc • (it — it') 
|fe • (it — it')| 2 + z/ 2 



z x ik' z x ik 



x- 



(u/, - o/jVcu-t =F a; ± W 



|fc? |fc |, t v- + -,-.. — -j (50) 

where we have used Eq. (1411) . It should be noted that the obtained diffusion fluxes can 
be divided into two parts, namely the diffusion tensor D(r,t) proportional to Vu± and the 
drift velocity V(r,t) proportional to u±. 



r ± = -D • Vw± ± Vu±, 



(51) 



(27 

7r5(k ■ (u — it')) 



ik ■ (u — it') 
\k ■ (u — it') I 2 + v 2 
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JS X ulxi Z X %r\j , . . . 

x , _ - (u;^ — w_) 



life 



/|2 



|fc| 2 

dk 



V = "/*'/ (2^/(2^ eXpW * + *>■<•—■'» 



(2?r) 2 y (2vr) 

vT(5(fc • (it — It')) — 



ik ■ (u — u') 



k ■ (u — u')\ 2 + v 2 



z x ik! z x ik 



\k 



/|2 



Ifcl 



w 



(52) 



(53) 



V. SPACE- AVERAGED COLLISION TERM 

Equation ( 150]) includes the oscillatory term exp(z(fc + fc') ■ (r — r')). To reveal the char- 
acteristics of the obtained collision term, we need to calculate the space average of the 
diffusion fluxes to drop the high-frequency component. Space average is calculated over 
the small rectangular area A(r) with sides both 2L located at r. The space average of the 
diffusion flux T± is defined by 



<r ± ) s = r s± (r) = -i-/ dr"T ± {r"). 
|A(r)| Jk{r) 



(54) 



We assume that the macroscopic variables such as u and u may be constant inside A(r) 
and only the term exp(i(fc + k') ■ (r — r')) should be space-averaged: 

(eMi(k + k')-(r-r'))) s 
= — — [ dx" f dy"exp(i(k + k') ■ r") exp(-i(k + k') ■ r') 

(2Ly j-L J-L 



7T 



-) 5(fc + fc')exp(-z(fc + fc')V) 



7T 



-) S(k + k') 



where r" = (x",y"). Thus, space-averaged diffusion fluxes are given by 



(55) 



r.±(r) 



- n (i)7^/^*-(«-^ 



z x fc z x fc 

|fc| 2 Ifcl 2 



(2tt)< 



(56) 



In Eqs. (|56|) . we omit the imaginary part as the collision term consists of only the real part. 
Further integration over fc in Eqs. (156]) can be performed. The integral concerning fc is as 



13 



follows: 

/dM(M«-«')) ( * > y . (57) 

Dividing A; into the parallel and the perpendicular components and inserting them into Eq. 

(JSZD, 

k = k\\n\\ + k±fi±, 
u — u' 



n» 



we obtain 



\u — u \ 
fix — z x nil, (58) 



(z x k)(z x fe) 



/ dkS(k ■ (u - u 1 )) 
J |fc 4 | 

dk\\ / dfcj_5(fc|||u — u'\) 

[z x (fcynn + A;j_nj_)][z x (ft||n|| + k±n±)] 



\k 2 + k 2 ± \ 2 

dk±- tttt(z x k±h ± )(z x fc±n±) 

| it — u'\ k\ 

J 1 u — u' 1 u — u' 

dk i — 



\u 



u'\ k 2 \u — u'\ k 2 , In 



2 



1 



(it — it') (it — it') 

i 7n z ~ / 

\u — u'\ 6 V k\ 

(u-u')(u-u') 2 

It It ft-min 



oc 



(59) 



where parameter fc min is introduced to regularize a singularity and determined by the largest 
wave length that does not exceed a system size, namely k m \ n = 27r/R where R is a charac- 
teristic system size determined by an initial distribution of the vortices. 

Finally, we obtain the following formulae for the diffusion and the drift terms. 

r,±(r) = -D, • Vw± ± V s oo ±} (60) 

s ~ (27r) 3 UJ k mm J \u-uf ' lbij 

Vs -{2<Kf\L) k mi J dr \u-uf • [bZ) 



It is worth stressing that the term V s u + has the opposite sign to the term V s 0J-. It may 
provide a mechanism for the "charge separation" which is usually seen in the equilibrium 
distribution for systems with positive and negative vortices J8|, [l8j. The notations r s± are 
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Q I 


dr' 


(u- 


- u 


')(«- 


«') 




RJ 




u 


-u'| 3 






K 


— U) 


-)V 


'jj - 


- (w+- 


- bJ_ 


.)V'u/ 



defined in the same manner as Eqs. (15 ip . In Eqs. (|6T|) and (162]) . two unknown parameters 
L and /c m i n remain. We assume that L = gi? where g is a size factor (j < 1). If we set 
g = 1/47T, total diffusion flux is rewritten as 

r s (r) = r s+ + r s _ 



(63) 



A. Diffusion flux in local and global equilibrium states 

At first, let us examine if the diffusion fluxes ( 1631) disappears in a local equilibrium state. 
We rewrite Eq. (163]) into a symbolic form: 

T s (r) = -^Jdr'j[u,^u',ij'] (64) 

where 7 is a functional of u, ip, u/ and ip'. Consider a state where temperature is locally 
uniform in each small region in the system. We call this state the local equilibrium state in 
which the local equilibrium condition is satisfied: 

wioq± = w ±exp(=F/3f2?/>ieq). (65) 

Inserting Eqs. ( 165]) into 7 in Eq. (164"]) . we find that 

_ pUleq ~ < q ) 

I '13 

I u lcq — W leql 

x (t*^ - u[ cq ) ■ [(u[ cq+ - uj' leq _) V(w lcq+ + w lcq _) - (w leq+ - w to q_)V'(wi B q+ + uj' leq _)} 
= -(3Q(u[ cq+ - u[ cq _)(u icq+ - w lcq _) 



x 



(ttteq - n; cq ) 

I '13 

l**leq — n leql 



>leq - U[ eq ) ■ (Wleq - V^Lq) 



= (66) 

where U\ cq = —z x V^i eq is used. As U\ cq — u[ cq is perpendicular to V^ieq — VVLq; 7 is equal 
to zero and this result indicates that the detailed balance is achieved. When the system 
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reaches a global thermal equilibrium state with uniform /3 |8| 

u eq ± = u ± exp(=F/3fiVeq), (67) 

15 



we obtain 



= -«< + - <_)(V% - V^cq + V^eq)- 

As (lieq — mL) • (V'^' — V^ eq ) = 0, the drift term in Eq. (1621) is rewritten as 



V s 



s,eq 



-/3fiD s , cq • V$ 



cq v if-"eq 



which is a counterpart of the Einstein relation lCfl . 



(68) 



(69) 



B. Energy-conservative property of collision term 



Time derivative of the total mean field energy E is given by 

dE 1 r t ,„. ,. (duj' ,dcu" 

— = - dr dr'Gir - r) -^-u + u/— 
dt 2 J J v ' \dt dt , 

a i duj 



where 



E = - / dri/xj 

= - fdrfdr'G(r-r')co'uj. 



Inserting the space-averaged equation of motion 

du 



Of 



+ V • (uu) 



v ■ r„ 



into Eq. (1701) . we obtain 
dE 



— = jdr^(-V-(uu)-V-T s ) 

/ drVip -uuj + drVip ■ T s 

f drVip ■ T s 

u — it') (it — u') 



— I dr J dr'Vi) 



(70) 



(71) 



(72) 



\u — u 



'13 



(u/ -u/jVw- (w + -w_)V / w / . (73) 



By permuting the dummy variables r and r' in Eq. ( 1731) and taking the half-sum of the 
resulting expressions, we obtain 

— = -— / dr / dr'(Vip - Vty 



df 



2R 



\u — u 



'|3 



x(u — it')- (a/ — wl)Vu; — (u + — w_)W 



0. 



(74) 
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We conclude that the obtained diffusion fluxes conserve the total mean field energy. 



C. H theorem 



The entropy function S is defined by using the H function: 

S = -k B H, (75) 

= — / drco + In uj + — w_ In |w_ | — 2N In Vt + const. (76) 

The time derivative of the H function is given by 

dH 1 r , du+ „ _ N duj. 



dt n 



_y dr _t (lnw+ + 1) (h |„.| + 1, 



— / druu + • V In w + — uw_ • V In |w_ | 
+— / drT s+ ■ V In w+ - T s _ • V In |u;_ 



— / dri* ■ Vw + — w ■ Vw_ 



n 

1 /• 

• V In w-l — r,_ • V In w_ 



= - fdrT s+ - Vlnw+-r s _- Vln|w_|. (77) 

For simplicity, we introduce the following notations. 

Vlnoj + = f + , Vln \oj~\ = v , 

u-u' = U. (78) 

Inserting Eq. (1631) into Eq. (1771) . we obtain 

lt = ~RJ J + ' |£/T ' [(W + ~ } + " U+ ] 



1 /• /" J7J7 

- dr dr'v^ ■ — -= ■ [(u' + - w^)Vw_ + w_W 
ft 7 j \U \ 



R 

-4 /dr /drS ] 



RJ J |C/ 3 | 
u) + u)' + {v + ■ UU • (v + — v' + ) 



+U-Lj'_(V- ■ UU ■ (v- — v'_ 



—uj+uj'_(v+ ■ UU ■ (v A 



J -\"+ 



v' 



uj_uj' + (v_ ■ UU ■ (t>_ + v' + )\ . (79) 
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By permuting the dummy variables r and r' in Eq. (!79|) and taking the half-sum of the 
resulting expressions, we obtain 



dt 2RJ J \U 3 \ 

x[uW+\(v + -v' + )-U\ 2 

+oj_oj'_\(v_-v'_)-U\ 2 

-u+u'_\(y+ - v'_) ■ U\ 2 

-uj_uj' + \(v_ - v' + ) ■ U\ 2 } <0. (80) 

The integrand of Eq. (1HD1 is positive or equal to zero, and dH/dt is negative or equal to 
zero. It is concluded that the entropy function (JTSJ) is the monotonically increasing function. 

VI. DISCUSSION 

We have demonstrated the simple and explicit formula of the Fokker-Planck type collision 
term for nonaxisymmetric, double-species point vortex system without the collective effect. 

Again, it should be noted that the drift term plays an important role in the inverse 
cascade process in 2D turbulence. In Eq. (160 p . the sign of the drift term changes following 
the sign of the vorticity, while the sign of the diffusion term is always negative. This implies 
that the drift term provides the "charge separation" of the vortices, which is commonly 
observed in equilibrium states at negative temperature. 

If the sign of (3 is negative, du eq /dip eq > 0, as 

duj P 



-"eq 



-f3Q(uj eq+ -w cq _) (81) 



#eq 

where 

^eq = ^cq+ + ^cq- (82) 

and the relations fIBTl) are used. The point vortex system easily approaches a quasi- stationary 
state near the thermal equilibrium state by a violent relaxation, which is purely collisionless 
and driven by mean field effects. The state is characterized by 

V ■ (uw) = (83) 

or equivalently 

lo = uj(iIj). (84) 
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In this state, we may expect that 



dip 



>0 



(85) 



almost everywhere in the system. If the above postulate is justified, the following two 
conclusions are drawn which elucidate the role of the drift term at negative temperature. 

Energy conservation is achieved by the energy dissipation process due to the diffusion 
term and the energy production process due to the drift term. We divide the expression 
( 1731) into two parts, namely the term which corresponds to the diffusion term and the one 
to the drift term. 



dE 
~dt 

dE 

~dl 

dE 

~dt 



dE 
~dt 
Q 
~R 



+ 



dE 
~dt 







v 



dr / drVip 



[u — u')(u — u') 



\u — u 



'1 3 



(u' + -u'_)Vu 



fi /" , [ , ,„ , (u-u')(u- u') 
— dr dr'Vip ■ - 



\u — u 



'13 



Wj 



If the vorticity u is a function of the stream function ip (see Eq. 
are rewritten as: 



dE 
~dt 

dE 
~dl 



— I dr I dr 
R 



,\Vip ■ (it— u 



A 12 



(86) 

(87) 

w_)W (88) 

), Eqs. dHB and 

dbj 



\u — u 



/|3 



Wj 



u 



R 



dr dr ,W.(«-u>)\>, 



Thus, if cLu /dip > 0, it is concluded that 

dE 

~dt 



\u — u 
dE 



'13 



UJ- 



dip 
do/ 

] ~chj? 



>0 



(89) 
(90) 

(91) 



Similarly, the diffusion term increases the entropy, while the drift term decreases it. We 

divide the expression ([77]) into two parts (the entropy function S is used instead of the H 

function). 

dS dS 
dt dt 

U V 

Vw. • D • Vu;_ 



dS_ 
+ dt 



>0 



(92) 



dS 
~dt 



dS_ 
~dt 



ks f r Va;+ • D • Vwj 

— / dr 

R J LU+ 



k B 
R 



dr I dr 



/|Vo; + ■ (u — u 



'^ iVw. 



it — u 



f\\2 



UJ + 



UJ. 



(93) 



v 



kb 
R 

ks_ 
R 



drV ■ Vw 



, / , ir-, (u — u')(u — u') . . 
dr / dr'Vuj ■ : J --— ■ VV 



\u — u 



!\3 



(94) 
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Equation (|93|) indicates that 



dS 
~dt 



> 



(95) 



regardless of the sign of du/dip. On the other hand, if the assumptions (1841 and (1851) are 
valid, Eq. (194"]) is rewritten as 



~dt 



k_B_ 
R 



-— dr dr 



\u — u 



/|3 



rf0 e^' 



and we obtain 



d* 7 



<0. 



It is also concluded that in the thermal equilibrium state with /3 < 0, 



dS_ 
~dl 



dS_ 



<0. 



(96) 



(97) 



(98) 



This may indicates that to form vortex clumps, a background distribution outside the clumps 
is necessary to dump the entropy. This conclusion is supported by the nonneutral plasma 



experiments 
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22J and the numerical simulation [181 ]. 



There are several outstanding issues remaining. First, the final formulae flBTl) and fl62|) 
include unknown parameters k min and L. Second, the integrals in Eqs. (16T|) for D and (16"2j) 
for V contain the divergent integrand, although combined terms T s ± = — D s ■ Vu± ± V s u± 
are regularized. A more rigorous justification will be needed for fixing the above two issues. 
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